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PROCEDE ET P1SPOS1T1F DE CALIBRATION 
ET DE CORRECTION DE N1VEAUX PE GRIS 

DOMAINE TECHNIQUE GENERAL. 

La presente invention concerne un procede de calibration d'appareil 
apte a acquerir une sequence d'images radiographiques et de correction 
d'images d'un objet a observer. 

Plus precisement, elle concerne un procede de calibration d'appareil 
apte a acquerir une sequence d'images radiographiques et de correction 
d'images d'un objet a observer en soustrayant a chaque image de la 
sequence la variation intempestive d'un niveau de gris entre les images 
d'une meme sequence. 
ETAT DE L'ART. 

II existe de nombreux appareils aptes a acquerir une sequence 
d'images radiographiques. De tels appareils sont notamment utilises dan ; s le 

,'r. 

domaine medical. 

Un exemple d'appareil de radiographie est schematiquement 
represents a la figure 1. Sur cette figure, un appareil 1 comportanKune 
potence 5 (d'extension sensiblement verticale) comporte un dispositif 3 
formant emetteur d'un faisceau radiographique 4 en direction d'une tablette 
7 comportant un detecteur 2. Le detecteur 2 comporte classiquement une 
matrice de photodiodes. Le detecteur 2 recupere les rayons X apres la 
traversee d'un objet a observer place dans I'espace d'observation 6 situe 
dans le faisceau 4, entre d'une part i'emetteur 3 et d'autre part le detecteur 
2 du dispositif 1. Les extensions de I'emetteur 3 et du detecteur 2 peuvent 
par exemple etre horizontales et perpendiculaires a la potence 5, mais 
toutes les directions d'observation peuvent etre possibles, notamment du 
fait de la rotation possible de ('ensemble forme par I'emetteur 3 et le 
detecteur 2 autour d'une axe d'extension sensiblement horizontale. Le 
dispositif 1 comporte egalement des moyens de traitement 8 permettant 
d'acquerir et de traiter une sequence d'images d'un objet issues du 
detecteur 2. 
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La difference d'absorption du rayonnement X radiographique par les 
differentes parties de I'objet a observer permet d'obtenir des 
renseignements sur la composition de I'objet. En effet, une image formee 
sur le detecteur 2 comporte differents niveaux de gris, dont on peut deduire 
5 une information. Ainsi, si I'objet a observer est une partie d'un corps humain 
par exemple, les os apparaTtront clairement sur Timage acquise par le 
detecteur et se detacheront nettement de la partie formee par des muscles. 

De tels appareils ne presentent cependant pas totalement 
satisfaction. 

1 0 En effet, de fagon generate, lors de I'acquisition d'une sequence 

damages radiographiques, on constate toujours une variation du niveau de 
gris entre des images successives de la sequence. Ainsi, la figure 2 montre 
que dans une succession d'images acquises aux temps t, f t 2 , t 3 , par 
exemple, le detecteur mesure une augmentation de la variation L du niveau 

15 de gris. Ce phenomene est du a une remanence des informations 
radiographique d'une image a I'autre, ce qui fait varier le niveau de gris 
entre les images. Sur la courbe de la figure 2, on a represents la variation L 
du niveau de gris de I'image courante au temps tj par rapport a une 
premiere image acquise au temps 0. La figure 2 montre egalement que la 

20 variation L de niveau de gris dans la sequence depend de I'epaisseur et de 
la composition de I'objet que Ton observe. Ainsi, pour un objet observe 
ayant une premiere epaisseur, on obtient la courbe f,, tandis que pour un 
objet ayant une deuxieme epaisseur differente de la premiere epaisseur, on 
obtient une courbe f 2 differente de la premiere. 

25 La variation L du niveau de gris est due principalement au piegeage 

de charges dans les photodiodes du detecteur. La variation L peut 
egalement avoir de nombreuses causes. II peut s'agir notamment d'une 
montee en temperature des differents elements du dispositif. 

Dans tous les cas, la variation du niveau de gris d'une image a I'autre 

30 fausse les mesures acquises par Tappareil. 

En effet, dans le cas d'une grande remanence, c'est-a-dire de 
variations tres importantes, on peut degrader grandement la qualite et 
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['interpretation des images acquises. En effet, il peut y avoir apparition sur 
les images de « fantomes », c'est-a-dire de superpositions d'images 
acquises precedemment sur une nouvelle image d'un objet. Quand la 
remanence est plus faible et ne provoque pas I'apparition de fantomes, les 
5 mesures effectuees sont egalement faussees par la variation intempestive 
du niveau de gris d'une image a I'autre. Les variations peuvent etre du 
meme ordre de grandeur que la dynamique en niveaux de gris du signal 
que Ton cherche a detecter. 

Certains procedes permettent de s'affranchir de cette remanence 

10 dans des applications particulieres avec des appareils particuliers. 

En effet, dans certaines applications, il est possible de determiner la 
loi selon laquelle la remanence diminue dans le temps et retrancher cette 
remanence dans les images de la sequence. La figure 3 montre qu'entre 
deux acquisitions d'images effectuees aux temps t, et t 2 , la remanence 

15 diminue selon une loi exponentielle que Ton peut determiner. Pour cela, on 
effectue entre deux instants d'acquisition au moins une mesure de noir 
correspondant aux niveaux de gris delivres par le detecteur en I'absence 
d'exposition par les rayons X. La mesure de noir s'effectue au temps t m par 
exemple et permet de determiner la valeur de la remanence L m pour cet 

20 instant. Par consequent, grace a la connaissance des valeurs de la 
remanence aux temps t, par exemple et t m , on peut deduire la loi de 
decroissance de la remanence. De la loi de la decroissance de la 
remanence entre deux acquisitions, on peut deduire la valeur de la 
remanence au temps t 2 de I'image acquise suivante et ainsi corriger les 

25 images acquises lors d'une sequence. 

On comprend que ce procede de correction de sequences d'images 
suppose ('utilisation de dispositifs de mesure de noirs. Or, de tels dispositifs 
ne sont pas toujours disponibles ni accessibles sur les dispositifs 
radiographiques de I'etat de I'art. Par consequent, dans la plupart des cas, il 

30 est pratiquement impossible de pouvoir corriger la variation de niveau de 
gris dans une sequence d'images successives. 
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PRESENTATION DE L'INVENTION. 

L'invention propose de pallierces inconvenients. 

Un but de Tinvention est de proposer un precede et un dispositif 
permettant de corriger la remanence dans une sequence dlmages 
5 radiographiques. Le procede et le dispositif doivent pouvoir s'affranchir de 
rutilisation de dispositifs de mesure de noirs. 

Un autre but de l'invention est de proposer un procede et un 
dispositif permettant de calibrer un appareil apte a acquerir une sequence 
d'images radiographiques. La calibration est effectuee afin de pouvoir etre 
10 utilisee pour corriger les effets de la variation du niveau de gris dans une 
sequence d'images radiographiques d'un objet a observer. 

Un autre but de l'invention est de proposer un procede et un 
dispositif de calibration et de correction d'images radiographiques 
applicables a tous les phenomenes repetables de variations, e'est-a-dire 
15 des phenomenes de variations qui se repetent d'une mesure d'une 
sequence d'acquisition a I'autre lorsque Ton garde les parametres de 
d'acquisition identiques d'une mesure d'une sequence d'acquisition a 
I'autre. 

A cet effet, l'invention propose un procede de calibration d'appareil 
20 apte a acquerir une sequence d'images radiographiques et de correction 
d'images d'un objet a observer, caracterise en ce qu'H comporte les etapes 
selon lesquelles : 

- N pour chaque image d'une sequence acquise par I'appareil et pour une 
frequence donnee d'acquisition de la sequence, on calibre I'appareil en 
25 determinant la valeur de la variation d'une moyenne des niveaux de gris 
dans au moins une zone d'interet de I'image courante d'au moins un 
dispositif de calibration, la variation etant determinee par rapport au niveau 
de gris moyen de la premiere image de la sequence dans chaque zone 
d'interet ; 

30 - on reitere la determination de ladite variation pour une serie de 
sequences d'images acquises avec des dispositifs de calibration conduisant 
a des premieres images de niveaux de gris moyens differents d'une 
sequence a I'autre ; 
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- on corrige chaque image cTune sequence d'images de I'objet a observer 
comportant des zones d'observation ayant differents niveaux de gris en 
soustrayant de I'innage courante la variation d'un niveau de gris par rapport 
a la premiere image de I'objet, la soustraction etant fonction du niveau de 

5 gris considere de chaque zone d'observation. 

L'invention est avantageusement completee par les caracteristiques 
suivantes, prises seules ou en une quelconque de leur combinaison 
techniquement possible : 

- on approxime une representation graphique ayant respectivement pour 
1 0 ordonnee et pour abscisse : 

- le rapport ayant au numerateur la variation du niveau de gris moyen de 
rimage courante de la sequence du dispositif de calibration par rapport au 
niveau de gris moyen de la premiere image et au denominateur le niveau de 
gris moyen de la premiere image ; et r . 

15 - les differents niveaux de gris moyens de chaque premiere image ; 
par une fonction dont on connaTt les caracteristiques ; 

- la fonction est une fonction affine ; _ t: 

- le niveau de gris que Ton corrige dans chaque zone de chaque image de 
I'objet a observer est le niveau de gris moyen ; 

20 - le niveau de gris que Ton corrige dans chaque zone de chaque image de 
I'objet a observer est le niveau de gris median ; 

- on corrige le niveau de gris d'au moins un pixel choisi dans chaque image 
de I'objet a observer ; 

- la calibration a lieu avant I'acquisition de la sequence d'images de I'objet a 
25 observer ; 

- chaque valeur du niveau de gris moyen de la serie de sequences de la 
calibration est donnee par I'observation d'au moins une plaque de 
calibration d'epaisseur determinee constituant chaque dispositif de 
calibration, son epaisseur changeant d'une serie a I'autre ; 

30 - on determine le niveau de gris moyen sur plusieurs zones d'interet 
simultanement lors de la calibration ; 
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- la soustraction lors de I'etape de correction depend d'une part du niveau 
de gris dans chaque zone d'observation et d'autre part de la position de la 
zone d'observation par rapport a chaque zone d'interet ; 

- la valeur soustraite a chaque image de la sequence damages de I'objet est 
5 fonction d'une part de la position de la zone d'observation et d'autre part 

d'une fonction de gain spatial determinee ; 

- la calibration a lieu pendant ('acquisition de la sequence d'images de 
I'objet a observer ; 

- chaque dispositif de calibration est place dans un champ d'acquisition de 
10 Tappareil comprenant egalement I'objet a observer ; 

- chaque dispositif de calibration comporte au moins deux zones d'interet 
ayant un niveau de gris moyen different d'une zone a I'autre pour chaque 
image ; 

- la valeur soustraite a chaque image de la sequence d'images de I'objet est 
15 fonction d'une part de la position de la zone d'observation et d'autre part 

d'une fonction de gain spatial de I'appareil ; 

- on combine les mesures de calibration d'au moins deux acquisitions de 
sequences successives ; et 

- on combine les mesures de calibration d'au moins deux acquisitions de 
20 sequences successives ayant des dispositifs de calibration differents. 

L'invention concerne egalement un appareil apte a acquerir une 
sequence d'images radiographiques mettant en oeuvre un precede selon 
l'invention. 

PRESENTATION DES FIGURES. 
25 D'autres caracteristiques, buts et avantages de l'invention 

ressortiront de la description qui suit, qui est purement illustrative et non 
limitative, et qui doit etre lue en regard des dessins annexes sur lesquels : 

- La figure 1, deja commentee, represente schematiquement un dispositif 
apte a acquerir une sequence d'images radiographiques selon I'etat de 

30 Tart ; 

- La figure 2, deja commentee, represente schematiquement revolution de 
la remanence en fonction du temps dans une sequence d'images 
radiographiques ; 
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- La figure 3, deja commentee, represente schernatiquement une methode 
de correction de la remanence selon I'etat de la technique utilisant une 
mesure de noirs ; 

- La figure 4, represente schernatiquement une premiere approche d'un 
5 precede possible de calibration et de correction de sequences d'images 

selon Tinvention, selon laquelle la calibration a lieu avant ['acquisition de la 
sequence d'images de I'objet a observer ; 

- La figure 5, represente schernatiquement une variante adaptative 
spatialement du precede selon la figure 4 ; 

10 - La figure 6 represente schernatiquement une deuxieme approche d'un 
procede selon Tinvention, selon laquelle la calibration a lieu en meme temps 
que I'acquisition de la sequence d'images de I'objet a observer ; 

- La figure 7 represente schernatiquement une variante adaptative 
spatialement du procede selon la figure 6 ; . 

15 - La figure 8 represente schernatiquement une coupe longitudinale d'un 
dispositif de calibration utilise dans la deuxieme approche ; et 

- La figure 9 represente schernatiquement les differentes etapes d'une 
variante possible de la deuxieme approche ou Ton combine des donnees de 
calibration d'au moins deux sequences successives ; et . 

20 - La figure 10 represente schernatiquement les differentes etapes d'une 
variante possible de la deuxieme approche ou Ton combine des donnees de 
calibration d'au moins deux sequences successives ayant des dispositifs de 
calibration differents. 

DESCRIPTION DETAILLEE DE L'INVENTION. 
25 L'invention propose un procede de calibration d'appareil apte a 

acquerir une sequence d'images radiographiques et de correction d'images 

d'un objet a observer permettant de corriger les effets indesirables des 

variations de niveaux de gris dans une sequence d'images radiographiques. 
Un procede selon l'invention peut etre considere selon deux 
30 approches. Selon une premiere approche possible, I'etape de calibration est 

effectuee avant I'acquisition de la sequence d'images de I'objet a observer. 

La deuxieme approche permet d'effectuer la calibration en meme temps 

que I'acquisition des images de I'objet a observer. 
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Dans toutes les figures, les elements ou etapes similaires sont 
references par des numeros identiques. 
Premiere approche. 

Voyons tout d'abord un procede possible mettant en ceuvre la 
5 premiere approche. 

Les principales etapes d'un procede possible utilisant la premiere 
approche sont representees schematiquement a la figure 4. Le procede 
represents a la figure 4 comporte deux grandes etapes principales. La 
premiere etape est I'etape de calibration. Elle se compose principalement 
0 de trois etapes references par (a), (b), et (c) sur la figure 4. La deuxieme 
grande etape du procede est la correction des images de I'objet a observer. 
Elle se compose principalement de trois etapes referencees par (d), (e) et 
(f). 

Decrivons tout d'abord I'etape de calibration. 

3 Lors d'une etape (a), on acquiert grace a I'appareil radiographique 

une sequence d'images referencees par #1, #2 #N. L'acquisition de la 

sequence est effectuee en observant un dispositif de calibration 40. Le 
dispositif de calibration 40 est positionne dans la zone d'observation au- 
dessus du detecteur de I'appareil et recouvre la majeure partie de la surface 

) du detecteur. 

Cette sequence est acquise pour une frequence d'acquisition 
donnee. Par exemple, la frequence d'acquisition peut etre d'une image 
toutes les 30 secondes, ou d'une image toutes les 60 secondes. Les 
images peuvent egalement etre acquises a intervalle de temps non regulier. 
■ La frequence d'acquisition de la sequence de calibration est 
preferentiellement la meme que la frequence d'acquisition que I'on souhaite 
utiliser pour l'acquisition de la sequence d'images de I'objet a observer. 

La figure 4(a) represente ainsi schematiquement une premiere 
sequence d'images d'au moins un dispositif de calibration 40 ayant une 
epaisseur donnee. Chaque dispositif 40 de calibration comporte au moins 
une plaque dont I'epaisseur est comprise entre 1 a 8 cm. Chaque plaque 
presente des caracteristiques d'absorption du rayonnement X emis par le 
dispositif sensiblement egales a I'objet que I'on veut observer dans la suite. 
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Ainsi, les plaques du type BR 12 connus de I'homme du metier sont 
preferees pour des utilisations en mammographie, puisqu'elles presentent 
les memes caracteristiques d'attenuation que les tissus glandulaires, 
comme le sein. Bien entendu, pourd'autres applications, d'autres materiaux 

5 peuvent etre utilises, comme par exemple la lucite ou le plexiglas. 

Lors de I'acquisition, on determine, pour chaque image, un niveau de 
gris moyen sur une zone d'interet choisie homogene et referencee par 41 
sur la figure 4. La dimension de la zone d'interet 41 est typiquement de 100 
pixels x 100 pixels. Pour ameliorer la precision des mesures du niveau de 

10 gris moyen dans la zone 41, on repete plusieurs fois I'acquisition des 
images du dispositif de calibration. On repete ainsi typiquement quatre a 
cinq fois chaque sequence. On moyenne ensuite les resultats. 

On effectue une serie d'acquisitions de sequences avec des 
dispositifs de calibration 40 ayant des epaisseurs differentes. On obtient 

15 ainsi des niveaux de gris moyens sur la zone d'interet 41 differents entre 
chaque sequence de la serie. 

On rappelle que le niveau de gris moyen evolue en fonction du 
temps, c'est-a-dire qu'il varie entre les images successives d'une meme 
sequence. On releve ainsi les niveaux de gris moyens dans la zone d'interet 

20 41 de chaque image de toutes les sequences de la serie. 

Ensuite, pour chaque nieme image C (n) d'une sequence acquise par 
I'appareil, et pour la zone d'interet 41 , on determine, grace a des moyens de 
traitement compris dans I'appareil, la valeurde la variation du niveau de gris 
moyen entre I'image courante et la premiere image. On note L (n) cette 

25 difference. Elle est egale a : 

Jin) (-,) 

On passe ensuite a I'etape referencee par (b) sur la figure 4. Dans 
I'etape (b), on determine pour chaque nieme image la difference relative 
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On trace ensuite le graphe de en fonction de C (1) . On 

constate alors que cette fonction peut etre approximee par une fonction 
affine. La difference estimee (l in) ) est alors egale a : 

(l m ) =C^x [a (/,) .C^ + /? (n) ] (2)ou a (n) et J3 M sont les coefficients 

5 de la regression lineaire des courbes tracees lors de I'etape (b) et calcules 
lors d'une etape (c). 

L'etape (c) montre que Ton determine, grace aux moyens de 
traitement compris dans I'appareil, les coefficients de regression a 0) et 
/? (0 correspondant respectivement aux coefficients directeurs et a 

10 I'ordonnee a I'origine de chaque droite de I'image i. On stocke ces 
coefficients dans des moyens memoire de I'appareil compris dans les 
moyens 8 de traitement ou dans des moyens situes hors de I'appareil. Cette 
derniere etape (c) termine la calibration. 

Chaque droite depend d'une part du rang de la nieme image dans la 

15 sequence, comme I'indique I'exposant (n), et d'autre part de la frequence 
d'acquisition de la sequence, Autrement dit, et comme on peut le constater 
sur la figure 4(b), pour une frequence d'acquisition identique, les courbes de 
I'image #2 et de I'image #3 sont differentes. Cela signifie egalement que 
pour une frequence d'acquisition differente, les courbes des images #2 sont 

20 differentes. 

On trace ainsi Tensemble des courbes correspondant a toutes les 
images de la sequence de calibration. 

On va maintenant decrire I'etape de correction des images de I'objet 
a observer. 

25 La figure 4(d) montre schematiquement que Ton acquiert une 

sequence d'images d'un objet a observer 42, la sequence comportant N 
images. Sans correction, on observe une variation du niveau de gris 
indesirable entre les images successives de la sequence. 

On va done utiliser les donnees de calibration afin de corriger ces 

30 variations indesirables. 
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Pour appliquer une correction au pixel de coordonnees (ij) d'une 
nieme image Y {n) non corrigee, on mesure le niveau de gris moyen Y$ wJ) 
dans une zone d'observation R(iJ) centree sur le point (ij) de la premiere 
image, par exemple la zone referencee par 43 sur la figure 4(d). Chaque 
5 zone d'observation R(i 9 j) a une taille typique de I'ordre de 25 pixels x 25 
pixels. La taille de chaque zone d'observation correspond sensiblement a la 
plus petite taille des objets que Ton veut observer, par exemple une tumeur 
dans le cas duplications de I'invention a la mammographie. 

On applique ensuite a I'image Y in) le precede de correction en 
10 soustrayant de I'image courante la variation d'un niveau de gris par rapport 
a la premiere image de I'objet. 

En effet, la figure 4(e) montre que grace a la determination du niveau 

de gris moyen Yff ijy lors de I'etape (d) et grace aux donnees de calibration, 
on peut remonter a la valeur de la variation qui est une fonctiori de 

15 7j!> ~ . II suffit alors de retrancher cette valeur a la valeur du niveau de gris 

courante. ? 

On effectue ainsi une mesure du niveau de gris moyen dans 
plusieurs zones d'observation, par exemple egalement dans la zone 
d'observation referencee par 44 sur la figure 4(d). Le niveau de gris moyen 
20 de la zone d'observation 44 peut etre different du niveau de gris moyen de 
la zone 43. 

Mathematiquement, on obtient une nieme image corrigee Z\ n ) par 
la formule : 

25 En appliquant une telle correction sur Tensemble des images de la 

sequence, on obtient une sequence dans laquelle toutes les variations de 
niveaux de gris moyens ont ete supprimees, ou tout au moins reduites. 
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Dans les formules qui precedent, on remarque que Ton a utilise la 
mesure du niveau moyen de gris autour du pixel (ij). Plusieurs variantes 
sont possibles en fonction du niveau de gris que Ton considere. 
Premiere variante. 

5 On peut egalement determiner la valeur du niveau median de gris 

dans chaque zone d'observation. On ne considere done plus le niveau de 
gris moyen. On rappelle que la valeur mediane d'une serie est la valeur 
situee au milieu de la serie de valeurs placees dans un ordre croissant ou 
decroissant. 

10 ^utilisation de la mediane permet d'eviter I'impact des valeurs du 

niveau de gris mesurees aberrantes sur la valeur du niveau de gris prise en 
compte dans Petape de correction. En effet, les valeurs extremes de la serie 
n'ont aucune influence sur le calcul de la valeur de la mediane. De telles 
valeurs aberrantes peuvent etre mesurees dans une zone de I'objet 42 

15 ayant un fort gradient d'epaisseur et aux alentours d'un changement brutal 
de I'epaisseur de I'objet a observer. 

On considere le niveau median de gris au lieu du niveau moyen de 
gris pour les objets presentant de tels gradients d'epaisseur. Les etapes de 
calibration et les formules (1 ), (2) restent identiques. Seule la formule (3) est 

20 legerement modifiee, puisqu'on I'applique au niveau de gris median et non 
plus au niveau de gris moyen. 

Le calcul de la mediane augmente generalement legerement le 
temps de traitement des informations effectuees dans les moyens de 
traitement de I'appareil. 

25 Deuxieme variante. 

Selon une deuxieme variante possible, le niveau de gris moyen peut 
egalement etre remplace par la valeur Y™ du niveau de gris du pixel (ij). 
On ne moyenne done plus la valeur sur une zone d'interet R{iJ) (43 et/ou 
44 par exemple) et Ton ne calcule plus la valeur mediane. 

30 De la meme fagon que pour la premiere variante, les etapes du 

procede restent identiques, et seule la formule (3) est legerement modifiee 
puisqu'on I'applique a la valeur Y$ du niveau de gris au pixel (ij). 
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La seconde variante permet, comme la premiere variante, d'avoir 
une bonne estimation de la variation du niveau de gris en des pixels situes 
pres d'une zone dans laquelle I'epaisseur de Tobjet a observer varie 
rapidement. Elle permet, contrairement a la premiere variante, une 
5 diminution du temps de traitement. Cependant, la determination de la valeur 
de la variation du niveau de gris est moins precise, puisqu'il y a 
augmentation des effets du bruit quantique. Les effets du bruit quantique 
etaient significativement reduits par le moyennage ou le caicul de la 
mediane. 

10 Dans tous les developpements qui precedent, on suppose que la 

variation du niveau de gris est homogene sur toute la surface du detecteur 
de I'appareil. On ne determine en effet les coefficients a {n) et j3 in) qu'en 
une seule zone d'interet 41 . 

APPROCHE ADAPTATIVE. 

15 On va maintenant decrire une approche adaptative de la premiere 

approche. L'approche adaptative prend en compte la non homogeneite des 
variations de niveau de gris sur la surface du detecteur de I'appareiL' En 
effet, il apparait que les detecteurs ne sont pas parfaits, et qu'il existe»:une 
disparite des variations de niveau de gris selon la position sur le detecteur. 

20 Les principals etapes d'une telle approche adaptative sont 

representees schematiquement a la figure 5. 

La methode adaptative consiste a prendre en compte la non 
homogeneite des variations en effectuant une calibration sur plusieurs 
zones d'interet 41. Les zones d'interet sont reparties de fagon reguliere sur 

25 la surface du dispositif de calibration et couvrent au maximum la surface du 
detecteur. 

Ainsi, lors d'une etape (a), on acquiert une sequence d'image d'un 
dispositif de calibration 40 dispose sur la surface du detecteur. Comme 
precedemment, on effectue acquisition d'une serie de sequences, afin de 
30 faire varier I'epaisseur du dispositif de calibration 40. Ensuite, pour chaque 
image C (/,) , on calcule la remanence L { £ par la formule : 




(similaire a la formule (1) deja mentionnee) 
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ou chaque R c ( k J\ (kj)<=Q represente une zone d'interet centree sur 
le pixel de coordonnees (k,l) referencee par 41, et ou Q. est I'ensemble des 
pixels sur lesquels sont centrees les zones d'interet referencees par 41 sur 
la figure 5. 

5 On calcule egalement z£ } /C (l) que \> on trace en f 0 nction de . 

On obtient ainsi lors de I'etape (b), non plus un graphe d'un 
ensemble de droites correspondant aux differentes images, mais un 
ensemble de graphes, correspondant chacun a une zone d'interet 41 . On a 
ainsi autant de graphes que de zones d'interets 41 sur la surface du 

0 dispositif de calibration 40. On a reference par 41 1 , 412,..., 41 N I'ensemble 
des graphes correspondants a I'ensemble des mesures dans les N zones 
d'interet 41. 

On approxime lors de I'etape (c) les courbes obtenues par des 
fonctions affines dont on determine, pour chaque graphe, les coefficients 
5 a™ et fi™ . Les coefficients a£> et fi™ dependent desormais de la zone 
d'interet 41 pour laquelle ils ont ete calcules. 

Lors de I'etape (d), on acquiert une sequence damages d'un objet 42 
a observer. 

On determine des zones d'observation 43 et 44 pour lesquelles on 
veut effectuer une correction. On repere alors la place des zones 
d'observation 43 et 44 par rapport aux differentes zones d'interet 41 . 

On prend alors, pour determiner la variation de niveau de gris de 
chaque zone d'observation, les valeurs de a£° et qui ont ete 

determinees pour la zone d'interet 41 correspondante. C'est ce qui est 
effectue lors de I'etape (e). 

On applique ensuite la correction grace a une formule similaire a la 
formule (3). Le niveau de gris corrige Z$ devient : 

z u = Wj + A ( ; } )J (similaire a la formule (3)) 

ou Yff est le niveau de gris avant la correction et a$ , fifties 

coefficients determines pour chaque zone d'interet et applique a la zone 
d'observation centree sur le pixel (ij). 
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On obtient ainsi a I'etape (f) une image Zf"J corrigee. 

Variantes possibles. 

Plusieurs variantes de I'approche adaptative sont possibles. 
Une premiere variante possible consiste a interpoler les coefficients 
a w e t p\ n ) de zones d'observation situees en dehors de zones d'interet. 

On peut dont ainsi avoir des coefficients aj") et fiff pour n'importe quelle 

zone d'observation. 

D'autres variantes sont egalement possibles et utilisent les valeurs 
medianes de niveau de gris ou les valeurs du niveau de gris en un pixel au 
lieu de considerer le niveau de gris moyen dans une zone d'observation. 

La premiere approche considere une calibration avant I'acquisition de la 
sequence d'images de I'objet a observer. Une deuxieme approche possible 
effectue la calibration en meme temps que 1'acquisition de la sequence 
15 d'images de I'objet. 
Deuxieme approche. 

Nous allons maintenant decrire une seconde approche d'unprocede 
selon I'invention, dans lequel les etapes de calibration ont lieu en meme 
temps que I'acquisition des images de I'objet a observer. 
20 Les principales etapes de la deuxieme approche sont representees 

schematiquement a la figure 6. 

On constate que sur la figure 6(a), le dispositif de calibration 40 est 
place dans un champ d'acquisition de I'appareil, et ce pendant I'acquisition 
de la sequence d'images de I'objet 42 a observer. 
25 On evite ainsi de devoir effectuer une serie d'acquisition uniquement 

pour la calibration. 

On a vu cependant qu'il etait necessaire dans la premiere approche 
d'effectuer une serie d'acquisitions avec des epaisseurs de dispositif de 
calibration 40 differentes. Cette variation d'epaisseur est necessaire au 
30 trace des courbes de I'etape (b). 

Dans la deuxieme approche, le dispositif de calibration 40 comporte 
au moins deux zones d'interet, ayant des niveaux de gris moyens differents 



5 



10 
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d'une zone a ('autre, pour chaque image. C'est ce que Ton a represents 
schematiquement a ia figure 8. En effet, sur cette figure, on constate que ie 
dispositif de calibration 40 comporte une premiere zone d'interet 41 
comportant des proprietes d'absorption des rayons X differentes d'une 
5 deuxieme zone 45. La difference de niveau d'absorption peut etre due a 
une difference d'epaisseur du dispositif 40 au niveau des deux zones 41 et 
45, et/ou a une difference de materiau des deux zones. Avantageusement, 
le dispositif comporte plus de deux zones. On comprend que plus le nombre 
de zones ayant des proprietes differentes est important, plus le trace lors de 

10 I'etape (b) des courbes de calibration aura un nombre de points significatifs 
important. Plus le trace des courbes en (b) est precis, plus leur 
approximation en (c) est precise egalement. 

Ainsi, lors d'une meme acquisition d'une sequence d'images du 
dispositif de calibration 40, on a plusieurs points de reference permettant 

15 d'effectuer le trace selon la figure 6(b) similaire a celui de la figure 4(b) de la 

premiere approche. On calcule pour chaque image Y {n) la remanence par la 
formule : 

L w = y(«) _ptij ( S i m ilaire a la formule (1) deja mentionnee) 
On trace ensuite les courbes de ces valeurs relatives par rapport au 
20 niveau de gris moyen dans la premiere image en fonction du niveau de gris 
dans la premiere image. II s'agit de I'etape (b). 

On approxime en (c) les differentes courbes des variations du niveau 
de gris en fonction du niveau de gris de la premiere image par des fonctions 
affines, et on en calcule les coefficients de la fonction representative. 
25 Lors de I'etape (e) on utilise la calibration des etapes (a), (b) et (c). 

La calibration permet de calculer la variation, pour chaque image de I'objet 
42 de la sequence. On obtient ainsi en (f) une sequence d'images 
corrigees. 

On comprend alors les avantages de cette approche. 
30 Premierement, on obtient une correction des variations du niveau de 

gris plus precise. En effet, la remanence est determinee directement a partir 
des images de la sequence de Tobjet a observer. On est done sur qu'il n'y a 
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pas de differences de comportement de Tappareil entre les sequences de 
calibration et Tacquisition des images de Tobjet a observer. 

Deuxiemement, ii n'est pas necessaire d'effectuer toute une serie de 
mesures afin d'obtenir les courbes de calibration. La calibration est 
5 effectuee directement, en meme temps que Tacquisition des images de 
Tobjet a observer. Ce deuxieme avantage represente un gain de temps 
considerable pour un utilisateur de I'appareiL 
Variantes de la deuxieme approche. 

On comprend qu'il n'est pas facile d'appliquer une approche 
10 adaptative a la deuxieme approche. En effet, du fait de la presence de 
Tobjet a observer, on ne peut effectuer une calibration sur Tensemble du 
champ d'acquisition. Ainsi, selon une premiere variante de la deuxieme 
approche, on utilise un modele spatial afin de tenir compte des eventuelles 
disparites de variations de niveau de gris en fonction de la position sur le 
15 detecteur. Les principales etapes d'une telle variante sont representees 
schematiquement a la figure 7. 

Les etapes (a), (b), (c) et (e) restent identiques a celle de la figure 6. 
On ajoute cependant entre Tetape (e) et Tetape (f) une etape (e') permettant 
Tapplication d'un modele spatial des disparites des variations. 
20 Ainsi, pour appliquer une correction a un pixel (ij) de Timage Y {n) , le 

niveau de gris moyen Y™ n est mesure dans une zone d'observation, 
referencee par 44 par exemple, centree autour du point (ij). 

Apres la calibration, le niveau Z/" } de gris du pixel (ij) devient 

25 (similaire a la formule (3)) 

avec Y { { f le niveau de gris avant la calibration, a in) et p w les 

coefficients de Tapproximation affine de la variation du niveau de gris et X i} 
un facteur de gain mesure ou modelise pour le pixel (ij). 

Pour obtenir un facteur de gain X u mesure, on peut par exemple 
30 utiliser une calibration pour plusieurs zones d'interet comme dans 
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I'approche adaptative decrite pour la premiere approche. Les differentes 
valeurs de X i j seront ainsi entrees dans les moyens de traitement de 
I'appareil et appliquees lors de I'etape de correction. 

Pour modeliser le facteur de gain X Lj , on peut utiliser I'inverse du 

facteur de gain de I'appareil. Le facteur de gain K jfi de I'appareil compense 
la non uniformite d'eclairement du detecteur par I'emetteur et la non 
uniformite de la reponse des photodiodes du detecteur. 

On observe done, sur le detecteur, une image corrigee en gain de la 
forme : 

Kjj (Im + remanence) ou Im est I'image de I'objet non corrigee en 
gain et non homogene. La partie (K i(j Im) est une image homogene. 

Ainsi, la remanence non homogene que I'on observe sur le detecteur 
est en fait : 

(Kjj remanence) 

La non homogeneite de la remanence observee est due au facteur 
de gain Kg. Ainsi, pour retrouver une remanence homogene, on applique au 
terme (Ky remanence) le facteur k i} tel que 



C'est la remanence homogene obtenue apres cette multiplication que 
Ton doit soustraire a I'image courante. 

Selon une deuxieme variante possible de la deuxieme approche, on 
compense le fait que peu de valeurs pour le trace des courbes de 
calibration peuvent etre acquises. En effet, comme on acquiert les donnees 
de la calibration en meme temps que les donnees d'observation de I'image, 
on n'a qu'une valeur de la variation par image pour un niveau de gris donne. 

Afin d'augmenter le nornbre de points de mesure pour ameliorer la 
precision du calcul des coefficients de regression a i0 et fi {0 de 
I'approximation, il est possible de combiner les informations de plusieurs 
sequences successives d'acquisition d'objets a observer differents. La 
figure 9 represente schematiquement les differentes etapes de cette 
variante. On apergoit sur la figure 9 qu'au moins deux sequences damages 
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d'objet 42' et 42" sont acquises lors d'etapes referencees par (ai) et (a 2 ). 
Les deux sequences sont referencees par 91 et 92, et sont acquises toutes 
les deux dans les memes conditions d'acquisition, mais successivement. 
Par exemple, la sequence 91 est acquise avant la sequence 92. La 
5 sequence 91 a notamment la meme frequence d'acquisition que la 
sequence 92. Lors de I'etape (b), les mesures recoltees sur le dispositif 40 
de calibration dans chacune des deux sequences 91 et 92 sont combinees 
afin d'augmenter le nombre de points sur lesquels on calcule des 
coefficients de regression lors de I'etape (c). Sur la figure 9, le dispositif 40 

10 comporte trois zones d'interet referencees par 41, 45 et 46. La figure 9(b) 
represente les differents points issus de chaque sequence, pour chaque 
image pour les zones 41, 45 et 46. On double ainsi le nombre de points si 
on combine deux sequences. La precision des approximations est done 
augmentee. Les deux sequences 91 et 92 sont ensuite corrigees lors des 

15 etapes (eO et (e 2 ) respectivement. La correction utilise les coefficients de 
regression calcules lors de I'etape (c). On obtient alors en (fi) et en (f 2 )'des 
images corrigees. 

La variante de la figure 9 est bien entendu generalisable a autant de 
sequences qu'on le desire. Pour chaque image, on a ainsi un nuage de 

20 points pour chaque niveau de gris. 

De plus, la variante de la figure 9 peut etre iterative, e'est-a-dire que 
pour chaque nouvelle sequence acquise, on peut recalculer des nouveaux 
coefficients de regression. 

Une autre variante de la deuxieme approche compense le fait que 

25 peu de points de fonctionnement sont disponibles pour le trace des courbes 
de calibration. En effet, du fait de la presence de I'objet a observer dans le 
champ d'acquisition du detecteur, on ne peut pas disposer dans le champ 
des dispositifs de calibration de grande taille. Les dispositifs ne possedent 
par consequent pas de nombreuses zones de caracteristiques d'attenuation 

30 differentes. La presente variante utilise des dispositifs de calibration 
differents d r une sequence a I'autre, afin d'augmenter le nombre de points de 
fonctionnement. Les principals etapes d'une telle variante sont 
representees schematiquement sur la figure 10. On apergoit sur la figure 10 
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qu'au moins deux sequences d'images d'objet 42' et 42" sont acquises lors 
d'etapes referencees par (ai) et (a 2 ). Les deux sequences sont referencees 
par 91 et 92, et sont acquises toutes les deux dans les memes conditions 
d'acquisition, successivement. La sequence 91 a notamment la meme 
frequence d'acquisition que la sequence 92. Les dispositifs de calibration 
40' et 40" utilises dans les deux sequences sont au contraire differents 
d'une sequence a I'autre. Par exemple, le dispositif 40' possedent deux 
zones d'interet 41 et 45, tandis que le dispositif 40" en possede au moins 
une, preferentiellement deux, referencees par 47 et 48. Les zones 47 et 48 
possedent des caracteristiques d'attenuation des rayons X differentes de 
celle des zones 41 et 45. Lors de I'etape (b), les mesures recoltees sur les 
dispositifs 40' et 40" de calibration dans chacune des deux sequences 91 et 
92 sont combinees afin d'augmenter le nombre de points de fonctionnement 
sur lesquels on calcule des coefficients de regression lors de I'etape (c). Les 
points de mesure des zones 47 et 48 lors de ('acquisition de la sequence 92 
completent les points des zones 41 et 45 acquis pendant I'acquisition de al 
sequence 91. On dispose par exemple les deux dispositifs 40' et 40" 
altemativement sur le detecteur lors d'acquisitions successives de deux 
objets differents 42' et 42". 

Dans les variantes des figures 9 et 10, on peut bien entendu utiliser 
un modele spatial. 

Bien entendu egalement, d'autres variantes de la deuxieme approche 
sont egalement possibles, et utilisent les valeurs medianes de niveau de 
gris ou les valeurs du niveau de gris en un pixel, au lieu de considerer le 
niveau de gris moyen dans une zone d 'observation. 

Generalisations. 

L'introduction du facteur de gain X u pendant I'etape de correction 
peut etre appliquee a la premiere approche egalement. 

Dans les developpements qui precedent, la variation estimee relative 
(Z- (n) )/C (l) dans la nieme image peut etre approximee une fonction affine du 

niveau de gris moyen dans la premiere image C (l) . Bien entendu, dans les 
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developpements qui precedent, on peut approximer la representation 
graphique visible a la figure 4(b) par exemple par d'autres fonctions dont on 
peut determiner ou dont on connaTt les caracteristiques. La fonction 
d'approximation choisie depend du phenomene repetable que Ton observe, 
5 ainsi que de la precision que Ton souhaite obtenir lors de la calibration. La 
fonction d'approximation peut ainsi par exemple etre une fonction 

polynomiale comportant des puissances de C (l) ou une fonction 
exponentielle. 

Le procede et I'appareil des developpements qui precedent 
10 s'applique avantageusement & des applications a la mammographie et a la 
detection de tumeurs cancereuses. 
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REVENDICATIONS. 

1 Procede de calibration d'appareil apte a acquerir une sequence 
d'images radiographiques et de correction d'images d'un objet a 
observer, caracterise en ce qu'il comporte les etapes selon 
lesquelles : 

pour chaque image d'une sequence acquise par I'appareil et pour une 
frequence donnee d'acquisition de la sequence, on calibre I'appareil en 
determinant la valeur de la variation d'une moyenne des niveaux de gris 
dans au moins une zone d'interet de I'image courante d'au moins un 
dispositif de calibration, la variation etant determinee par rapport au 
niveau de gris moyen de la premiere image de la sequence dans chaque 
zone d'interet ; 

on reitere la determination de ladite variation pour une serie de 
sequences d'images acquises avec des dispositifs de calibration 
conduisant a des premieres images de niveaux de gris moyens 
differents d'une sequence a I'autre ; 

on corrige chaque image d'une sequence d'images de I'objet a observer 
comportant des zones d'observation ayant differents niveaux de gris en 
soustrayant de I'image courante la variation d'un niveau de gris par 
rapport a la premiere image de I'objet, la soustraction etant fonction du 
niveau de gris considere de chaque zone d'observation. 

2 Procede selon la revendication 1, caracterise en ce qu'on approxime 
une representation graphique ayant respectivement pour ordonnee 
et pour abscisse : 

le rapport ayant au numerateur la variation du niveau de gris moyen de 
I'image courante de la sequence du dispositif de calibration par rapport 
au niveau de gris moyen de la premiere image et au denominateur le 
niveau de gris moyen de la premiere image ; et 
les differents niveaux de gris moyens de chaque premiere image ; 
par une fonction dont on connait les caracteristiques. 
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3 Procede selon la revendication 2, caracterise en ce que la fonction 
est tine fonction affine. 

4 Procede selon Tune des revendications 1 a 3, caracterise en ce que 
le niveau de gris que Ton corrige dans chaque zone de chaque 
image de I'objet a observer est le niveau de gris moyen. 

5 Procede selon Tune des revendications 1 a 3 f caracterise en ce que 
le niveau de gris que Ton corrige dans chaque zone de chaque 
image de I'objet a observer est le niveau de gris median. 

6 Procede selon Tune des revendications 1 a 3, caracterise en ce 
qu'on corrige le niveau de gris d'au moins un pixel choisi dans 
chaque image de I'objet a observer. 

7 Procede selon I'une des revendications 1 a 6, caracterise en ce que 
la calibration a lieu avant ('acquisition de la sequence d'images de 
I'objet a observer. 

8 Procede selon la revendication 7, caracterise en ce que chaque 
valeur du niveau de gris moyen de la serie de sequences de la 
calibration est donnee par ('observation d'au moins une plaque de 
calibration d'epaisseur determinee constituant chaque dispositif de 
calibration, son epaisseur changeant d'une serie a i'autre. 

9 Proc6de selon I'une des revendications 7 ou 8, caracterise en ce 
qu'on determine le niveau de gris moyen sur plusieurs zones 
d'interet simultanement lors de la calibration. 

10 Procede selon la revendication 9, caracterise en ce que la 
soustraction lors de i'etape de correction depend d'une part du 
niveau de gris dans chaque zone d'observation et d'autre part de la 
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position de la zone d'observation par rapport a chaque zone 
d'interet. 

1 1 Procede selon Tune des revendications 7 ou 8, caracterise en ce 
que ia vaieur soustraite a chaque image de la sequence damages de 
I'objet est fonction d'une part de la position de la zone d'observation 
et d'autre part d'une fonction de gain spatial determinee. 

12 Procede selon Tune des revendications 1 a 6, caracterise en ce que 
la calibration a lieu pendant ('acquisition de la sequence d'images de 
I'objet a observer. 

13 Procede selon la revendication 12, caracterise en ce que chaque 
dispositif de calibration est place dans un champ ^'acquisition de 
I'appareil comprenant egalement I'objet a observer. 

14 Procede selon la revendication 13, caracterise en ce que chaque 
dispositif de calibration comporte au moins deux zones d'interet 
ayant un niveau de gris moyen different d'une zone a I'autre pour 
chaque image. 

15 Procede selon Tune des revendications 12 a 14, caracterise en ce 
que la vaieur soustraite a chaque image de !a sequence d'images de 
I'objet est fonction d'une part de la position de la zone d'observation 
et d'autre part d'une fonction de gain spatial de I'appareil. 

16 Procede selon Tune des revendications 12 a 15, caracterise en ce 
qu'on combine les mesures de calibration d'au moins deux 
acquisitions de sequences successives. 



17 Procede selon I'une des revendications 12 a 16, caracterise en ce 
qu'on combine les mesures de calibration d'au moins deux 
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acquisitions de sequences successives ayant des dispositifs de 
calibration differents. 

18 Appareil (1) apte a acquerir une sequence d'images radiographiques 
d'un objet (42) a observer, caracterise en ce qu'il comporte des 
moyens (8) aptes a : 

pour chaque image d'une sequence acquise par I'appareil et pour une 
frequence donnee d'acquisition de la sequence, calibrer I'appareil 
comportant des moyens aptes a determiner la valeur de la variation 
d'une moyenne des niveaux de gris dans au moins une zone d'interet 
choisie de I'image courante d'au moins un dispositif de calibration (40), 
et determiner ladite variation par rapport au niveau de gris moyen de la 
premiere image de la sequence dans chaque zone d'interet ; 
reiterer la determination de ladite variation pour une serie de sequences 
damages de dispositifs de calibration acquises avec des dispositifs de 
calibration conduisant a des premieres images de niveaux de gris 
moyens differents d'une sequence a I'autre ; 

corriger chaque image d'une sequence d'images de I'objet (42) a 
observer comportant des zones d'observation ayant differents niveaux 
de gris et a soustraire de I'image courante la variation d'un niveau de 
gris par rapport a la premiere image de I'objet, et a effectuer la 
soustraction en fonction du niveau de gris considere de chaque zone 
d'observation. 

19 Appareil (1) selon la revendication 18, caracterise en ce qu'il 
comporte des moyens (8) aptes a approximer une representation 
graphique ayant respectivement pour ordonnee et pour abscisse : 

le rapport ayant au numerateur la variation du niveau de gris moyen de 
I'image courante de la sequence du dispositif de calibration par rapport 
au niveau de gris moyen de la premiere image et au denominateur le 
niveau de gris moyen de la premiere image ; et 
les differents niveaux de gris moyens de chaque premiere image ; 
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par une fonction dont un utilisateur de I'appareil (1) peut determiner 
les caracteristiques. 
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